var c h k;
parameters A beta theta delta G J kss;
beta = 0.98;
theta = 0.36;
A = 0.5;
delta = 0.1;
G = (1/(beta*theta)-(1-delta)/theta)^(1/(1-theta));
J = 1/(beta*theta)-(1-delta)/theta - delta;
kss = (1-theta)*(delta+J)/(G*(A*J+(1-theta)*(delta+J)));

model;
A*c/(1-h) = (1-theta)*k^theta*h^(-theta);
c(+1)/c = beta*(theta*k^(theta-1)*h^(1-theta)+1-delta);
k = (1-delta)*k(-1)+k(-1)^theta*h^(1-theta)-c;
end;

initval;
k = kss;
h = G*kss;
c = J*kss;
end;

steady;
resid;